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ABSTRACT 

Mult i- wavelength observations of the black hole X-ray binary XTE J1118+480 
have offered abundant spectral and timing information about the source, and have 
thus provided serious challenges to theoretical models. We propose a coupled 
accretion-jet model to interpret the observations. We model the accretion flow 
as an outer standard thin accretion disk truncated at a transition radius by 
an inner hot accretion flow. The accretion flow accounts for the observed UV 
and X-ray emission, but it substantially under-predicts the radio and infrared 
fluxes, even after we allow for nonthermal electrons in the hot flow. We attribute 
the latter components to a jet. We model the jet emission by means of the 
internal shock scenario which is widely employed for gamma-ray bursts. In our 
accretion-jet model of XTE Jll 18+480, the jet dominates the radio and infrared 
emission, the thin disk dominates the UV emission, and the hot flow produces 
most of the X-ray emission. The optical emission has contributions from all three 
components: jet, thin disk, and hot flow. The model qualitatively accounts for 
timing features, such as the intriguing positive and negative time lags between the 
optical and X-ray emission, and the wavelength-dependent variability amplitude. 



Subject headings: accretion, accretion disks — black hole physics — ISM: jets 
and outflows — stars: individual (XTE J1118 + 480) — X-rays: stars 
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1. Introduction 

Strong evidence now exists for black hole primaries in 15 X-ray novae (also known as soft 
X-ray transients; McClintock & Remillard 2004). One such source — XTE Jl 118+480 — was 
discovered with the All-Sky Monitor aboard the Rossi X-Ray Timing Explorer (RXTE) on 
2000 March 29 (Remillard et al. 2000). Subsequent optical observations led to a measurement 
of the mass function, f{M) — 6.00 ± O.36M0, which represents a lower limit on the mass 
of the compact primary and thus makes the source a secure black hole candidate (BHC; 
McClintock et al. 2001a; Wagner et al. 2001). XTE J1118+480 is one of the best observed 
BHCs. It lies at an unusually high Galactic latitude (-1-62°), close to the "Lockman Hole" 
region. The foreground absorption is extremely low (with Nh ~ 0.7 — 1.3 x 10^°cm~^; 
Hynes et al. 2000; McClintock et al. 2001b), which allowed the detection of the source by 
the EUVE satellite (Hynes et al. 2000). Simultaneous (or near-simultaneous) observations 
were conducted, on multiple occasions, at radio, infrared, optical, UV, EUV, and X-ray 
wavelengths, with state-of-the-art instruments (Hynes et al. 2000; McChntock et al. 2001b; 
Prontera et al. 2001; Chaty et al. 2003; McClintock et al. 2003). 

For clarity, we briefly summarize the main observational results here. These include 
two aspects — spectral and timing features. The most complete spectral energy distribution 
(SED) of XTE J1118+480 is shown in Figures 1 and 2. The radio data are from Fender et 
al. (2001) and the infrared to X-ray data from McChntock et al. (2001) (all the data are 
associated with "epoch 2", when the best simultaneous coverages were achieved; see Chaty 
et al. 2003 for a summary of all observations). The radio spectrum is well described by 
a power-law of the form F^, oc z/°'^. Such a spectrum is often thought to be typical of jet 
emission, although no jet has been directly imaged, down to a limit of < 65 D(kpc) AU 
(Fender et al. 2001), where D is the distance to the source. Note that we do not include in 
Figures 1 and 2 an observational data point at 350 GHz (Fender et al. 2001), because this 
measurement was not done simultaneously with the others. From IR to UV, the spectrum 
is flat, with the HST spectrum exhibiting emission lines. Also, a Balmer jump is seen in 
absorption ai u 10^^-^ Hz (Hynes et al. 2000), implying that thermal emission contributes 
substantially to the optical/UV band. The derived EUV spectrum depends sensitively on 
the assumed Nh, which is still not well constrained but probably lies in the range Nh = 
1.0 - 1.3 X lO^^cm"^ (McClintock et al. 2001b, 2004). We take this uncertainty into account 
by requiring the model to stay within the allowed range at EUV energies. McClintock et al. 
(2001b) fitted the X-ray spectrum with a broken power- law. Above ~ 2 keV they obtained 
a photon index of 1.78, but below ~ 2 keV the spectrum appeared to be relatively harder. 
However, cahbration issues were subsequently noted for the ACTS detectors used in the 
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Chandra observations^. This makes the spectrum uncertain at low energies. There is, in 
fact, independent evidence that the break at 2 keV may not be real. XTE J1118+480 was 
observed many times with BeppoSAX, but the X-ray spectra show no apparent deviation 
from a single power-law at low energies (Frontera et al. 2001). 

The main timing features include the following. 1) A quasi-periodic oscillation (QPO) 

feature was detected in the X-ray light curve, initially at a frequency u ~ 0.08Hz (Rcvnitsev, 
Sunyaev, & Borozdin 2000), and was subsequently found to evolve (Wood ct al. 2000). The 
QPO was also detected in the optical and UV bands at similar frequencies (Haswell et 
al. 2000; Yamaoka, Ueda & Dotani 2000). The fractional rms amplitude of the QPO is 
8 — 10% in the X-ray but only about 1% at UV wavelengths (Hynes et al. 2003, hereafter 
H03). The fact that the same QPO frequency is seen at optical, UV, and X-ray wavelengths 
indicates a common origin. 2) XTE J1118-I-480 also shows rapid aperiodic variability at most 
wavelengths. The variability amplitude is quite large both in the X-ray and IR bands but is 
small in the optical/UV band. 3) Correlation between emission at different wavelengths is 
apparent (H03). In particular, cross-correlation analysis has revealed some puzzling details 
in the correlation between the optical and X-ray emission (Kanbach et al. 2001; H03; Malzac 
et al. 2003). In general, the optical photons appear to lag the X-ray photons by 1 — 2 s (see 
H03, though with caveats). The lags are wavelength dependent; on average a longer delay 
is seen at longer wavelengths. On the other hand, the cross-correlation function (CCF) also 
shows a "precognition dip", i.e., the optical emission decreases about 2 — 5 seconds before the 
corresponding X-ray increase (Kanbach et al. 2001). At UV wavelengths the "dip" appears 
to be weaker and the lag becomes shorter, ~ 0.5 s (H03). These complicated positive and 
negative time lags between optical/UV and X-ray emission are not easy to understand. What 
is quite clear from the derived autocorrelation functions (ACFs) is that the optical/UV 
emission is not consistent with being due to the re-processing of X-ray photons by the 
accretion disk, as is often assumed, because the ACF at optical/UV wavelengths is narrower 
than that in X-rays (Kanbach et al. 2001; Spruit & Kanbach 2002; H03). 

Several models have been proposed to explain the observed spectral and temporal prop- 
erties of XTE J1118+480. Esin et al. (2001, hereafter EOl) explain the spectrum with an 
advection-dominated accretion flow (ADAF) model, based on the work of Narayan (1996) 
and Esin, McClintock, & Narayan (1997). They assume that the gas lost from the secondary 
initially forms a standard thin disk outside a transition radius rtr. At rtr, the cool disk is 
truncated and makes a transition to a hot accretion flow, described as an ADAF (Narayan 
& Yi 1994, 1995b; Narayan, Mahadevan & Quataert 1998). EOl satisfactorily explain the 
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X-ray, EUV and UV spectra of the source, but their model shghtly under-predicts the optical 
flux and significantly under-predicts the IR fluxes. They do not include radio measurements 
in their work, but it is quite clear that their model cannot account for the emission at radio 
wavelengths. 

In contrast, Markoff, Falcke, & Fender (2001) propose that the SED of XTE J1118+480 

is dominated by synchrotron radiation from a jet, although they also need a truncated 
accretion disk to explain the UV and EUV spectra. Inside the truncation radius, they 
assume that the accretion fiow becomes an ADAF-like accretion flow. However, unlike EOl, 
they ignore the radiation from the ADAF. 

No attempts have been made to explain the observed timing properties with either of 
the above models. Merloni, Di Matteo & Fabian (2000) consider both spectral and timing 
data in their work, but their magnetic flare model predicts that the disk emission should 
peak at about 0.2 keV, which is in disagreement with the EUVE and Chandra data. Also, 
the model implies almost no time lag between optical and X-ray photons, which seems to be 
at odds with the measurements. Recently, Malzac, Merloni & Fabian (2004) have proposed 
a time dependent, coupled disk-jet model for XTE J1118-I-480, which has some resemblance 
to the model we discuss in this paper. Whereas our model attempts to fit the spectral data 
(sec the following sections), Malzac et al. concentrate on understanding the timing features. 
As pointed out by them, due to the complexity of the time evolution of the accretion-jet 
system, detailed modeling is impossible. They thus adopt a phcnomcnological approach. 
They model the variability by assuming random fiuctuations of the output power from the 
disk and the jet, with the power being injected from a reservoir of stored magnetic field. 
By carefully choosing their parameters, they are able to reproduce almost all the observed 
timing features. These parameters can, in principle, constrain the dynamics and geometry of 
the accretion flow. One of their interesting results is that they can rule out models in which 
the energy budget is completely dominated by either the jet or the accretion flow; rather, 
they favor a model in which both components contribute. 

In the present paper, we describe a coupled accretion-jet model to simultaneously ac- 
count for both the spectral and timing properties of XTE J1118-I-480. We propose that the 
X-ray spectrum is produced mainly by the ADAF-like hot accretion flow, whereas the radia- 
tion at longer wavelengths comes from a jet (as in AGN). A similar idea has been suggested 
previously (e.g., Hynes et al. 2000; McClintock et al. 2001; Chaty et al. 2003). In § 2, we 
describe the model and discuss how it can explain the SED of XTE Jl 118+480. In § 3, we 
show that the observed temporal properties can also be accommodated qualitatively within 
the model. We conclude in §4 with a summary and discussion. We present in the Appendix 
technical details on calculating the jet emission. 
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2. Fitting the Spectrum 

2.1. Accretion flow 

The accretion component of our model is implemented in nearly the same manner as 
in EOl, i.e., the accretion flow consists of an inner ADAF and an outer thin disk. However, 
we have taken into account advances in our understanding of the ADAF during the past 
ten years. First, both numerical simulations (Stone, Pringle, & Begelman 1999; Hawley & 
Balbus 2002; Igumcnshchev et al. 2003) and analytical work (Narayan & Yi 1994, 1995a; 
Blandford & Begelman 1999; Narayan et al. 2000; Quataert & Gruzinov 2000) indicate that 
probably only a fraction of the gas that is available at large radius actually accretes onto 
the black hole. The rest of the gas is either ejected from the flow or is prevented from being 
accreted by convective motions. The details are likely to depend on the accretion rate. 

We note that the outflow (and convection) is ultimately the result of the accreting gas 
acquiring a positive Bernoulli parameter, as emphasized by Narayan & Yi (1994, 1995a). 
Further, the effect is strongest when the accretion rate is much below the threshold above 
which ADAF ceases to exist. Thus, accretion flows in highly under-luminous sources, like 
Sgr A* or quiescent X-ray binaries, are expected to have strong outflows. On the other hand, 
the Bernoulli parameter decreases with increasing radiative efficiency, and in fact becomes 
negative when the radiative cflicicncy is large enough. Therefore, for more luminous systems 
like XTE Jl 118+480 in outburst and other X-ray binaries in the low/hard state, which have 
relatively high accretion rates and radiate fairly efficiently, we expect outflows and convection 
to be less well-developed. In the present paper, we allow for this effect by adopting the 
following phenomenological prescription for the change in mass accretion rate as a function 
of radius. We assume that, in the hot flow, 

where 

s(r) = so/(r), (if 0</(r)<l), (2a) 

s{r) = 0, (if /(r) < 0). (2b) 

Here sq is a constant, which we set to sq = 0.3, as suggested by our previous modeling 
of the highly advection-dominated source Sgr A* (Yuan, Quataert & Narayan 2003). The 
parameter /(r) is the advection factor of the accretion flow, deflned as 

/W = 7- = ^^' (3) 

yvis yvis 

where ^adv; Qvis ^-nd Qie are the rates of energy advection, viscous heating, and Coulomb 
collision cooling for the ions, respectively. When the accretion rate is very low, as in the 
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case of Sgr A*, gvis S> Qie, so /(r) = 1 and s{r) — sq. In this case, from eq. (1) we have 
the usual form, M = M^ir / r^^Y^ , where Mq is the accretion rate at the transition radius rtr 
(or the outer boundary of the ADAF). We adopt Sq = 0.3, as in the case of Sgr A* because 
the physics of the outflow should be the same as long as /(r) = 1 even though the accretion 
rates (in Eddington units) can be quite different. We should note, however, that our results 
are not sensitive to the exact value of Sq. 

A negative value of / in eq. (2b) means that advcction plays a heating rather than a 
cooling role. In this case, the hot accretion flow is described by a luminous hot accretion flow 
(hereafter LHAF) model, which is a natural extension of an ADAF to higher accretion rates 
(Yuan 2001, 2003). From ADAF to LHAF, both M and the radiative efficiency increase 
continuously and smoothly. Yuan & Zdziarski (2004) argue that for luminous X-ray sources, 
such as the low/hard states of some BHCs and Seyfert 1 galaxies, the luminosity may be 
above the highest luminosity an ADAF can reach but could be accommodated by an LHAF. 
We allow for an LHAF in this work, because it is unclear at present which solution, ADAF or 
LHAF, applies to XTE J1118+480. We simply refer to both the ADAF and LHAF solutions 
as hot accretion flows. 

We calculate the global solution of the hot accretion flow, starting at rtr and integrating 
inward. The numerical details may be found in Yuan (2001). One main difference with 
EOl is that we solve the radiation hydrodynamics equations self-consistently, and thus we 
obtain the exact value of /(r) at each radius. In contrast, EOl used the approximation that 
/(r) has a constant average value at all radii. On the other hand, we treat Comptonization 
within a local approximation, whereas EOl computed the Comptonization globally using 
the method described in Narayan, Barret & McClintock (1997). The radiation processes 
we consider include bremsstrahlung, synchrotron emission, and the Comptonization of both 
synchrotron photons from the hot accretion flow and soft photons from the cool disk outside 
rtr. The emission from the outer cool disk is modeled as a multicolor blackbody spectrum. 
The effective temperature as a function of radius is determined by the viscous dissipation 
and the irradiation of the disk by the inner hot flow. 

Yuan & Zdziarski (2004) found that to explain the X-ray emission of most black hole 
X-ray binaries, a > 0.1 is required (see also Narayan 1996). We flx a and the magnetic 
parameter f3 (deflned as the ratio of the gas pressure to the sum of gas and magnetic pressure) 
at their "typical" values: a = 0.3, (3 = 0.9. We set 6 = 0.5, i.e., 50% of the viscous dissipation 
heats electrons directly. The exact value of 6 does not affect our results very much since the 
required M to model XTE Jll 18+480 in outburst is relatively high, so the main heating 
mechanism for electrons is energy transfer from ions via Coulomb coUisions. In this sense 
a, (5 and 5 are not free parameters, though we should emphasize that large uncertainties 



- 7- 



exist here. We set the mass of the black hole at M = 8Mq, the distance to the source at 
D = 1.8 kpc, and the binary inclination 9 = 70° (McClintock et al. 2001a; Wagner et al. 
2001). Following EOl, wc estimate the outer radius of the cool disk using Paczyhski's formula 
(Paczyhski 1971): rout = 3 x lOV,(lOM0/M)2/3, where r, = 2GM/c^ is the Schwarzschild 
radius of the black hole. The free parameters of the accretion flow arc the transition radius 
rtr, the accretion rate at the transition radius Mq, and an outer boundary condition — the 
temperature of the accretion flow at rtr (Yuan 1999). 

Figure 1 shows the spectral fitting results obtained with the accretion fiow model. The 
values of the parameters arc: Mq = O.OSMEdd, ^tr = 300r,5. The X-ray emission is produced 
by Comptonization in the hot fiow. The main seed photons are from synchrotron emission 
by the thermal electrons in the hot flow (as assumed in the original ADAF model of Narayan 
& Yi 1995b), as opposed to the blackbody emission of the thin disk. This is also consistent 
with the prediction of Wardzinski & Zdziarski (2000) given that XTE J1118+480 is not very 
luminous. For more luminous sources, the seed photons may be dominated by blackbody 
emission from the thin disk. The EUV and UV in the model are mostly from the outer thin 
disk. The flt is satisfactory, although the optical fluxes are slightly under-predicted. The 
fact that the UV/optical emission is dominated by the thin disk explains the presence of 
Balmer jump absorption and emission lines and reprocessing features in the data (§1). The 
IR and radio fluxes are signiflcantly under- predicted, however (ref. Fig. 2). Figure 3 shows 
the proflles of the advection factor /(r) and the fractional mass accretion rate M{r)/Mo as 
a function of radii. We see that /(r) is positive over much of the flow except near rtr- Since 
most of the radiation comes from the inner region where /(r) > 0, the solution is in the 
ADAF rather than LHAF regime, consistent with EOl. This is because the luminosity of 
XTE Jl 118+480 is not high. 

While our results are in general agreement with those of EOl, there are two noteworthy 
differences. First, our value of rtr(= 300r5) is signiflcantly larger than that of EOl (rtr = 
55rs). This discrepancy is mainly due to two reasons. First, EOl adopted a no-torque 
boundary condition at rtr while we apply this condition at the marginally stable orbit of 
the black hole. Second, in EOl the mass accretion rate of the thin disk follows M(r) = 
Mo(l — rtr/r) while we simply use M(r) = Mq. Both differences are related to the physics 
of the transition of the accretion flow at rtr, which is highly uncertain at present, so it is 
not clear which approach is more appropriate. As a comparison, rtr = 352r5 in Chaty et 
al. (2003) who fltted the EUV spectrum, while rtr = ITvs in Prontera et al. (2001; 2003) 
who fltted the iron line and reflection features. The second difference between our model 
and EOl is that the value of M in EOl (Mq = 0.02MEdd) is signiflcantly smaller than ours 
(Afo = O.OSMEdd)- This is primarily because (1) we include an outflow in our calculations 
so that the accretion rate in the inner region is smaller than that at rtr (see Fig. 3 where 
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M ~ 0.03MEdd near the black hole in our model, close to EOl's value); and (2) we use 

the pseudo-Newtonian potential of Paczynski & Wiita (1980), while EOl used the general 
relativistic solution of Popham & Gammic (1998) in calculating the radial velocity of the 
accretion flow. As shown by Narayan et al. (1998), the latter gives higher luminosity for the 
same accretion rate. 

To account for the under-prediction of the IR and radio fluxes, we first consider the effect 
of nonthermal electrons in the hot accretion flow. Since the inflowing gas is collisionless, 
processes such as MHD turbulence, reconnection, and weak shocks can accelerate electrons 
and generate a nonthermal tail at high energies in the electron distribution function. Yuan, 
Quataert & Narayan (2003) found that the radio spectrum of Sgr A*, which was under- 
predicted by a pure ADAF model with only thermal electrons, can be explained if roughly 1% 
of the electron energy is in nonthermal electrons. We tested this idea for XTE J1118-I-480. 
The dotted line in Figure 2 shows the (absorbed) synchrotron emission from nonthermal 
electrons. We see that there is a sharp cut-off below about 10^^ Hz, so that the emission 
from nonthermal electrons is unable to fit the radio and IR fluxes. This result is not sensitive 
to how much energy the nonthermal electrons have. In the case of Sgr A*, the emission from 
nonthermal electrons extends to much lower frequency and forms a power-law spectrum. 
The difference between Sgr A* and XTE Jl 1184-480 is that in the latter case the density 
is several orders of magnitude higher. Therefore, the magnetic field in XTE J1118-I-480 is 
much stronger and the lowest frequency that the power-law electrons emit is much higher. 
We conclude that the accretion fiow alone cannot account for the low-frequency spectrum of 
XTE Jl 118-1-480 at radio and IR wavelengths. Some other component, most hkely a jet, is 
required. 



2.2. Coupled Accretion-Jet Model 

Jets are thought to occur in the low/hard state of BHCs (see Fender 2004 for a review). 
There have been many papers on the emission of radio jets in active galactic nuclei (e.g., 
Blandford & Konigl 1979; GhiseUini, Maraschi, & Treves 1985; Falcke 1996). In the present 
paper, following Spada et al. (2001), we adopt the internal shock scenario widely used in 
interpreting gamma-ray burst (GRB) afterglows (e.g., Piran 1999). The details of the model 
of the jet radiation are described in Appendix A. Briefly, we assume that, near the black hole, 
a fraction of the accretion flow is transferred into the vertical direction to form a jet. Since 
the radial velocity of the accretion flow near the black hole is supersonic, a standing shock 
should occur at the bottom of the jet due to the bending. From the shock jump conditions, 
we calculate the properties of the postshock fiow, such as the electron temperature Tg. We 
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assume a constant Tg in the jet, which is clearly over-simplified, since adiabatic expansion 
will cause the electrons to cool. However, the assumption has very little effect on the results 
because the jet emission is dominated by the nonthermal electrons discussed below. We 
assume that the jet has a conical geometry with half opening angle 0, and that the bulk 
Lorenz factor of the jet Fj is independent of distance from the black hole. We further 
assume that internal shocks occur due to the collision of shells with different Fj. These 
shocks accelerate a fraction of the electrons into a power-law energy distribution with index 
p — 2.24 (e.g.. Kirk et al. 2000). The steady state energy distribution of the accelerated 
electrons is carefully determined since it is important for calculating the emitted spectrum. 
The effect of radiative cooling is considered in this process. Following the widely adopted 
approach in the study of GRBs, we specify the energy density of accelerated electrons and 
amplified magnetic field by two free parameters, and e^. We then calculate the radiative 
transfer by both thermal and power-law electrons in the jet, although we find that the latter 
plays a dominant role. Only synchrotron emission is considered since Compton scattering is 
not important in this case (sec also Markoff, Fender & Falcke 2001). 

The thin solid line in Figure 2 shows the emission of the jet. The parameters are: mass 
loss rate in the jet Mjct = 2.5 x 10^'^MEddi which is about 0.5% of the accretion rate in 
the accretion disk, = 0.1, — 0.06, = 0.02, bulk Lorenz factor of the jet Fj = 1.2, 
and length of the jet ~ 13 AU. The values of and eb are well within the typical range 
obtained in GRB afterglows (e.g., Panaitescu & Kumar 2001; 2002), and the length of the 
jet is consistent with the observed upper limit of 65£)(kpc) AU. The value of Fj is well within 
the range obtained by combining observations and numerical simulations: Fj < 1.67 (Gallo, 
Fender & Pooley 2003). We see from Figure 2 that the jet emission fits the low-frequency 
radiation very well. The IR fiux is dominated by the jet, while from optical to UV, the jet 
becomes less important. The contribution of the jet to EUV and X-rays is negligible. We 
should point out that the solution shown is not unique and that the jet parameters arc not 
as well constrained as those of the accretion fiow. However, the results are not very sensitive 
to the values of the jet parameters. 

It is interesting to check whether a pure thermal jet can also explain the data. We find 
that we can get an equally good fit to the spectrum if we adjust the geometry and T^^z) 
profile of the jet carefully. In this model, we only need a tiny fraction of the gas in the 
accretion flow, ~ 0.003%, to go into the jet. However, the required temperature is very high, 
Tg ~ 10^° K. In addition, the jet velocity has to be very low, ~ 100 km s~^; otherwise, the 
required magnetic fleld in the jet becomes unrealistically large. Such a low speed close to 
the black hole seems unphysical. 
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3. Interpreting the Timing Features 
3.1. QPOs 

Numerous models have been proposed to explain the QPO phenomenon in X-ray binaries 
(see review by van der Khs 2000). In some models, the QPO frequency is associated with the 
Keplerian frequency of the accretion flow at a special radius — the transition radius rtr in our 

case. For example, Giannios & Spruit (2004; see also RezzoUa et al. 2003) suggest that the 
QPO can be excited by the interaction of the inner hot accretion flow and outer thin disk. 
The QPOs then result from the basic p-mode oscillations of the inner hot accretion flow, 
with frequency roughly equal to the Keplerian frequency at rtr. The Keplerian frequency 
at rtr = 300rs is ~ 0.22 Hz, which is roughly consistent with the observed QPO frequency 
of ~ 0.1 Hz. Because the entire region of the hot flow oscillates collectively at the same 
frequency, and the emission from the hot flow contributes somewhat at both optical/UV and 
X-ray (see Fig. 1), the QPO should be observable at both optical/UV and X-ray wavelengths 
with the same frequency. Wood et al. (2000) find that the QPO frequency in XTE J1118-F480 
increases from 0.07 to 0.15 Hz during the outburst, while the 2-10 keV X-ray flux slowly 
rises and then decreases. Our calculations do not show such a non-monotonic relationship, 
so the evolution of the QPO remains a puzzle. We should emphasize that the non-monotonic 
change of the QPO frequency with the flux is not universal among BHCs. In fact, for most 
sources, the correlation seems to be monotonic (e.g., Cui et al. 1999). 



3.2. Variability amplitude 

The variability amplitude from the jet is expected to be large, both from internal shocks 
and from possible instabilities in the jet. The hot accretion flow is thermally marginally 
unstable, so any perturbations in it will survive and move inward, as shown by numerical 
simulations (Manmoto et al. 1996) and analytical work (Yuan 2003). However, the growth 
timescale of the perturbations is longer than the accretion timescale, so the hot accretion 
flow is not threatened by the instability. The simulations further show that the simulated 
flux variation can account for the observed substantial variability observed in BHCs. On 
the other hand, the intrinsic variability of emission from the thin disk should be very weak 
because the characteristic timescale is many hours even at rtr, i-e., much longer than the 
observed ~ seconds or minutes variability timescale (e.g., Kanbach et al. 2001). The only 
source of variability of the thin disk emission is due to the reprocessing of the variable X-ray 
radiation, but the contribution of this component is very weak. 

With the above knowledge, we can qualitatively understand variability amplitudes at 
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different wavelengths. Large variability in the IR and X-ray bands is natural because the IR 
emission is dominated by the jet and the X-ray emission by the hot flow. As the emission 
from the disk becomes more important in the optical and UV, the source varies less in these 
bands. The correlation between optical/UV and X-ray is easily understood because the 
hot accretion flow contributes in both bands. H03 find that the spectral energy distribution 
(SED) of the variable component of the emission is roughly a power law, which they argued as 
being consistent with optically thin synchrotron radiation. However, given the fact that the 
rms amplitudes were derived from light curves with the same time resolution, it is actually 
not straightforward to interpret the result, since the intrinsic variability timescale at different 
wavelengths should be quite different. Moreover, the physical origin of the variability is likely 
to be complicated (e.g., Malzac, Merloni, & Fabian 2004). We note that a power-law SED 
of the variability does not arise naturally in a pure jet model (e.g., Markoff, Falcke & Fender 
1999). For instance, if we assume that the variability is caused by fiuctuations in Mjct, such 
a model would predict a power-law index of 0.8, which is the same as the X-ray spectral 
index, while the measured index of the variability spectrum is ~ 0.59 (H03). 



3.3. Correlations between optical/UV and X-ray 

Suppose there is a perturbation due to an instantaneous increase of Mq. The X-ray 
flux will increase. The increase in M will propagate inward with the accretion flow, and 
eventually will lead to an increase in the mass loss rate and thus the optical/UV emission 
from the jet. This could explain why the optical/UV variability lags the X-ray variability. 
Quantitatively, we find that in our model the optical/UV emission from the jet comes mainly 
from regions at a distance of about d ~ GOOOr^ from the black hole. This corresponds to a 
propagation time of~d/c~1.2s, consistent with the measured ~ 1 — 2 s lag. The size of the 
optical emission region is ~ 2d(f) ~ 1200rs, where is the half opening angle of the jet. The 
corresponding light crossing time is 1200rs/c ~ 0.1s, consistent with the shortest variability 
timescale ~ 100 ms seen in the optical (e.g., Kanbach et al. 2001). Since the emission at 
longer wavelengths originates from regions farther away, the time lag should increase with 
increasing wavelength. 

As for the negative lag, we note that, for the parameters of our model (Fig. 1), an 
increase of M in the hot accretion fiow results in a decrease of the optical/UV fiux, as shown 
in Figure 4. The optical/UV emission from the hot accretion flow is mainly due to self- 
absorbed synchrotron emission, which depends on the proflles of Tg and optical depth r. 
For our model, an increase in M causes a decrease in the fiux. In our model, the optical 
emission comes from ~ SSr^, the UV from < lOr^, and the X-rays from ~ 7 — Sr^. So when 
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Mo increases, the optical flux wiU flrst decrease, then the UV wiU decrease, and flnaUy the 
X-ray flux will increase. This might be the origin of the negative lag of the optical/UV, as 
well as the negative correlation, and may also explain why the lag in the UV is shorter than 
in the optical. Since the emission from the hot accretion flow contributes less at shorter 
wavelengths in the optical/UV regime (see Fig. 1), we can also understand why the dip 
becomes weaker at shorter wavelengths. Since the IR flux from the hot accretion flow does 
not vary with varying Mq (see Fig. 4), we predict that such a negative lag should be absent 
between IR and X-ray. 

Quantitatively, however, wc arc not able to account for the magnitude of the negative 
lags. The viscous timescale at ~ SSr^ is ~ 0.1 s, which is more than 20 times smaller than 
the observed 2 — 5 s negative lag seen in the optical. This might be due to an approximation 
in the outer boundary condition we assume for the global solution. For technical reasons, we 
set the angular velocity of the flow at rtr to be substantially sub-Keplerian, il(rtr) ~ 0.5Qk, 
even though it should be super-Keplerian (Abramowicz, Igumenshchev, & Lasota 1998); 
otherwise, the viscous dissipation would be negative, and the solution would be unphysical 
(see also Manmoto, Mineshige, & Kusunose 1997). Since the centrifugal force is the dominant 
factor determining the radial velocity of the accretion flow, our approximation makes the 
radial velocity much larger than it should actually be and thus lead to a shorter viscous 
timescale. In addition, the viscosity parameter a may be smaller than the value we adopted, 
which will again result in a longer viscous timescale. 

Finally, we note that an increase of M in the cool thin disk will obviously result in an 
increase in the optical/UV emission. However, such an increase is unlikely to be seen in the 
cross-correlation analysis, since the accretion timescale in the thin disk is on the order of 
hours. 



4. Summciry and Discussion 

The observational data on XTE J 111 8+480 is almost unique among all current BHCs. 
The spectral and timing information impose very strong constraints on theoretical models 
and provide us with an opportunity to understand in detail the inflow/outflow processes 
around black holes. In this paper we explain how these observations can be understood in 
the context of a coupled accretion-jet model. In our model, the accretion flow is described 
as a geometrically thin cool disk outside a transition radius rtr and a geometrically-thick hot 
accretion flow inside rtr, as in the model of EOl. We adopt a phenomenological prescription 
for the magnitude of the mass outflow from the hot accretion flow (eqs. 1-3). The free 
parameters describing the accretion flow are the transition radius rtr, the mass accretion 
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rate at rtr, Mo, and the outer boundary condition at rtr- The spectrum due to the accretion 
flow alone is shown in Figure 1. The X-ray emission is dominated by Comptonization of 
synchrotron photons in the hot accretion flow, and both the EUV and UV arc dominated 
by the cool disk. The fit is quite satisfactory in these bands. The optical flux is slightly 
under-predicted, however, and the IR and radio spectra are signiflcantly under-predicted 
(Fig. 2). These results are very similar to those of EOl. 

Obviously, we require an additional component in the model to explain the IR and radio 
fluxes. We flrst consider the possibility of nonthermal electrons in the hot accretion flow, 
but flnd that this idea docs not work. Wc stress, however, that the failure does not mean 
that there arc no non-thermal electrons in hot accretion flows. Such electrons might, for 
instance, be responsible for the "hard tail" in the spectrum of Cyg X-1 in the low/hard state 
(McConnell et al. 2000). 

Having eliminated non-thermal electrons as an explanation for the low frequency emis- 
sion of XTE J1118-I-480, we argue that the radiation must originate in a jet. Assuming 
that a small fraction of the mass in the accretion flow is transferred to the jet, we calculate 
the jet emission using the internal shock scenario that is widely adopted in the study of 
GRB afterglows. The results of the accretion-jet model are shown in Figure 2. We find that 
the radiation from the jet can account for all of the radio and IR emission and part of the 
optical/UV emission. The required mass loss rate in the jets is about 0.5% of the accreted 
matter. 

The coupled accretion-jet model not only explains the spectrum, it also qualitatively 
explains many of the timing features observed in XTE J1118-I-480. These features include 
the frequency of QPO; the similarity of the QPO frequency in optical/UV/X-ray bands 
(§3.1); the dependence of the variability amplitude on wavelength (§3.2); and the positive 
and negative time lags between optical/UV and X-ray (§3.3). Quantitatively, however, we 
are not able to account for the magnitude of the negative time lag between X-ray and 
optical/UV (§3.3). 

It is interesting to examine the energetics of the accretion fiow and the jet in our model. 
The total accretion power is Pace = Mqc^ ~ 5 x 10^^ ergs~^ and the power lost in the outfiow 
is -Poutflow = -Pace — M{rs)c^ — 3.6 X lO^^ergs"^. The X-ray luminosity emitted by the hot 
accretion fiow is -Lx-ray ~ 2 x lO^^ergs"^, the optical/UV luminosity emitted by the thin 
disk is ~ 2 X lO^^ergs"^, and the jet power is Pjet = FjMjetC^ ~ 3.6 x lO^^ergs"^, which 
is ~ 2 times Lx-ray For comparison, Malzac et al. (2004) require P^et/L^-iay ~ 10 to 
reproduce the main timing features of XTE J1118+480, while Fender et al. (2001) estimate 
Pjct/-^x-ray ^ 0.2. The lumiuosity emitted by the jet in our model is Ljct ~ 2 x lO^^ergs"^, 
so the radiative efficiency of the jet is ~ 0.055, roughly consistent with the estimate of ~ 0.05 
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by Fender et al. (2001) but larger than the value of ~ 0.003 in Malzac et al. (2004). So 
there are differences in both the value of Pjet/Lx-ray ^-nd the efficiency of the jet between 
our model and that of Malzac et al. (2004). One reason for the discrepancy is that Malzac 
et al. assume the optical flux to be completely dominated by synchrotron emission from 
the jet, while our detailed modeling shows that the contribution from the accretion flow 
and the jet are comparable in the optical band (Fig. 2). Thus, more power from the jet is 
required in their model. In addition, the estimated value of Ljet in Malzac et al. (2004) is 
only 5 X lO^^ergs"^, which is ~ 4 times smaller than ours. This is because they integrate 
the jet emission from radio to optical, while in our model, the jet emission extends up to 
X-rays (ref. Fig. 2). 

Assuming Pj^c — -Poutflow ~ 1-4 x lO^^ergs"^ to be the accretion power in the inner 
region of the accretion flow from which most of the X-ray and jet power originate, we see 
that only i^x-ray/(-Pacc — -Poutflow) ~ 1% IS released through the X-ray emission and ~ 2% 
channeled into the jet, while most of the accretion power is stored in the accretion flow and 
advected into the black hole. In other words, XTE J1118+480 is radiatively quite inefficient, 
in agreement with the conclusion of Malzac et al. (2004). The small ratio of the jet power 
to the accretion power also justifies our approximation that the jet has very little effect on 
the global solution of the hot accretion flow. We should point out that some uncertainties 
exist in the above estimations concerning the jet since the jet parameters in our model are 
not as well constrained as the parameters of the accretion flow. 

Several other caveats also need to be mentioned. First, we adopt a pseudo-Newtonian 
potential rather than the exact general relativistic approach when we calculate the dynam- 
ics of the hot accretion flow. Secondly, we adopt a sub-Keplerian angular velocity at the 
transition radius whereas the rotation here should be super-Keplerian. The main effect of 
these two approximations is that the radial velocity in the hot flow is larger than it should 
actually be, and thus the density is smaller than the "correct value" . We believe that most 
of the effect is absorbed in the accretion rate parameter Mq. But the approximations do 
affect some quantitative result such as the time lag between optical/UV and X-ray. Thirdly, 
we have not explored fully the parameter space. The values of several parameters such as 
a, (3, and 6 are flxed in our calculations (to 0.3, 0.9 and 0.5, respectively). Investigating their 
effects in detail by surveying their entire parameter space would be very time-consuming 
and is beyond the scope of the paper. 

The philosophy of this paper is that the hard X-ray emission comes from the hot accre- 
tion flow via thermal Comptonization, and that the contribution from the jet is negligible 
in this band. This is different from the model of Markoff et al. (2001) in which synchrotron 
radiation from the jet dominates in X-rays. We note that many details of the X-ray ob- 
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servations of BHCs have been successfully explained with a hot accretion flow model (see 
the review by Zdziarski & Gierlihski 2004) and it remains an open question whether the 
jet model can do equally well. Poutancn & Zdziarski (2002) and Zdziarski et al. (2003) 
have pointed out some difficulties with the jet proposal. For example, the non-thermal syn- 
chrotron emission in this model cannot produce a sharp enough cut off at high energies, and 
the predicted spectrum is not as hard as the spectra observed in many BHCs. Also, the 
jet model should yield X-ray variability virtually independent of energy, which is in strong 
disagreement with the observational data. Finally, it is unclear if the model can explain the 
various timing features of XTE J1118-I-480 described in this paper. 

Of course for some black hole sources, the emission from the jet dominates over the 
accretion flow in the X-ray band. BL Lacs are a well-known class of objects where this 
situation is known to exist. In previous work we have discussed this possibility also for two 
other sources, Sgr A* and NGC 4258 (Yuan, Markoff, & Falcke 2002; Yuan et al. 2002). In 
the case of NGC 4258, the jet emission dominates the accretion flow because we require a 
significant fraction of the accretion fiow to be transfcrcd to the jet, Mjet/Mo ~ 10 — 25%, 
which is more than ~ 20 times higher than in XTE J 111 8+480. Such a high value perhaps 
implies that the black hole in NGC 4258 is very rapidly spinning. In the case of Sgr A*, 
the value of Mjet/Mo is similar to XTE J1118H-480, but the X-ray emission from the jet is 
comparable to the accretion flow. This is because the accretion rate (in Eddington units) 
in Sgr A* is much lower. The flux from the accretion flow, which comes from (multi-order 
scattering) Comptonization radiation, increases much faster with the accretion rate than that 
from the jet, which is from synchrotron and (one-order scattering) synchrotron-self-Compton 
emission. Therefore, the ratio of jet to disk flux increases with decreasing Eddington-scaled 
accretion rate. 

Recently a very interesting correlation between radio and X-ray fluxes has been dis- 
covered in GX 339-4. The correlation extends over more than three decades in X-ray flux 
(Corbel et al. 2003). Such a correlation likely exists in other BHCs and even in AGN (Gallo, 
Fender, & Pooley 2003; Merloni, Heinz & Di Mattco 2003; Falcke, Kording, & Markoff 
2004). The correlation is sometimes used as evidence for a jet origin for the X-ray emission 
of BHCs, e.g., Markoff et al. (2003). However, Heinz (2004; see also Merloni, Heinz & Di 
Matteo 2003) recently pointed out that if the electron energy spectrum is not too steep and 
if radiative losses are included, both of which are required by observations, the jet model 
cannot explain the radio — X-ray correlation. Merloni, Heinz & Di Matteo (2003) further 
showed that the X-ray emission is unlikely to be produced by radiatively efficient accretion 
(as in the sandwiched corona+disk geometry); rather, the accretion flow must be radiatively 
inefficient. Our preliminary investigations indicate that the radio — X-ray correlation can be 
explained in the context of our accretion-jet model (Yuan & Cui 2004, in preparation). 
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A. The Internal Shock Model for Jet Radiation 

We adopt the internal shock scenario to calculate the emission from the jet, similar 
to Spada et al. (2001). We are interested only in the time- averaged spectrum. Following 
Blandford & Konigl (1979), we assume the jet is in conical geometry, with semi angle of 
whose axis makes an angle 6 with the direction of observer. The jet has a constant velocity, 
characterized by a bulk Lorenz factor of Fj , and has constant plasma temperature. The mass 
loss rate in the jet is, 

Mjet = 7r^^(/)V(-2)^j (Al) 

The quantity p{z) is the mass density of the jet plasma at distance z from the black hole, 
measured in the jet-comoving frame. 

The main assumption in the internal shock scenario is that the central power engine 
produces energy which is channelled into jets in an intermittent way, thus faster shells will 
catch up with slower ones and internal shocks are formed in the jet. The minimum distance 
the shells propagate before collision occurs is zq ~ F^r^ (Piran 1999; Spada et al. 2001). 
Our results are not sensitive to its exact value. 

The bulk Lorenz factor of steady jets in BHCs is hkely only mildly relativistic (Fender 
2004), e.g., Fj < 1.67 from Gallo, Fender & Pooley 2003. In this case, for an adiabatic index 
of 4/3, the energy density of the internal shock is (Piran 1999), 

62 = '^2n2mp(? (A2) 

where 72 = ^^J (F? + l)/2 is the Lorenz factor of the formed internal shock, n2 — (472 + 3)ni is 
the post-shock number density with n\ is the preshock number density in the jet determined 
byeq. (Al). 

The shock will heat plasma in the jet, generate/amplify the magnetic field, and accel- 
erate a small fraction of electrons into relativistic energy. We assume that the fraction of 
accelerated electrons in the shock is and fix = 1%- Given the uncertainty in shock 
physics, as the usual approach, we introduce two dimensionless parameters, Cg and e^, which 
measure the fraction of the comoving internal energy of the internal shock stored in the 
accelerated electrons and magnetic field. Obviously, and are not independent. 

Assume that the injected electrons after the shock acceleration have a power-law distri- 
bution with index p, 

n^\{.le)d-ie = N^\{P - l)7e~^> 7min < 7e < 7max (A3) 

We set p — 2.24, according to the results of relativistic shock acceleration of Bednarz & 
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Ostrowski (1998) and Kirk et al. (2000). In this case (p > 2), we have 

^Pl = Cen^T^n ■ (A4) 
Now we calculate the value of 7min- We have, 

iVpl^eC^^3^7m7n = ^eUsh = ee(72 " 1)^2771^0^ (A5) 

where Ugh = (72 ~ l)n2mpC^ is the internal energy density of the internal shock. From the 
above equation and the definition of ^e, we can obtain 

7min = 72 - 1 7 T 0^6) 

p - 1 me 

The value of 7max is not important if we are not interested in the fitting the X-ray spectrum 
of XTE Jl 118+480 with jet emission. When radiative cooling of relativistic electrons is 
important, as in the present case of XTE J1118+480, the steady distribution of electrons 
is different from eq. (A3). Defining a "coohng Lorenz factor" 7c at which the radiative 
timescale trad is equal to the dynamical timescale tdyn at distance z in the jet, 

3 SnrrieC z , , 

then depending on the relative value of 7niin and 7c, there will be two cases for the steady 
distribution. When 7inin > 7c, we have, 

n^\{le)d-ie = N^\{P - I)7c7m7n7e"^^^7e, 7c < 7e < 7min, (A8a) 

npl(7e)c?7e = iVpi(p - l)7c7-(^+')ci7e, 7e > 7min. (A8b) 

When 7inin < 7c) we have, 

n^\{le)d'^e = Npi{p - l)j~Pdje, 7min < 7e < 7c, (A9a) 
nMdje = Nplip - l)7c7e"^^+'^C?7e, 7e > 7c- (A9b) 



The magnetic field generated/amplified by the shock is determined by, 

— = esUsh = €5(72 - l)n2mpC^ (AlO) 

OTT 



Since most of electrons may still be in thermal distribution, we need to consider their 
role in emitting and absorbing photons. To this purpose, we need to know their temperature. 
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One constraint comes from the following consideration. If the jet is formed at the innermost 
region of the accretion flow, within the sonic point at ~ lOrg, since the accretion flow is 
supersonic, when it is bended into the vertical direction to form the jet, a standing shock 
should occur. Note that the global solution of ADAF (e.g., Narayan, Kato & Honma 1997) 
does not find shocks. Our assumption of the bending shock is not in confiict with this result 
since jet was not considered in that calculation. On the other hand, shock is found in the 
general relativistic MHD numerical simulations of jet formation (e.g., Koide et al. 2000). 
^Prom the global solution of the accretion flow, we know the values of preshock quantities. 
Applying the shock jump conditions at the jet radius, we then be able to calculate the 
postshock quantities, including the electron temperature (see Yuan, Markoff, & Falcke 2002 
for details). Adiabatic expansion will cause the electrons to cool while the internal shocks 
in the jet will further heat the electrons. But for simplicity, we do not consider these effects, 
since we find the radiation from the power-law electrons dominate over that from thermal 
ones. 

Now we are ready to calculate the emission from the jet. The emissivity from each 
location in the jet is. 



where r is the optical depth along the line of sight in the jet, S^, = (jth + ipi)/(ctth + ctpi) 
is the source function, including the emission and absorption from both thermal (jth, ctth) 
and power-law {jpi, ctpi) electrons in the jet. We then integrate the emission from different 
distance in the jet to obtain the total emission. The relativistic effects is taken into account 
in the calculation. There is a remaining important point when we do the integration, that 
is, we should not integrate all of the volume of the jet. A "volume filling factor" /sh(< 1) 
should be introduced. The value of /sh is very uncertain. It obviously depends on the "spatial 
density" of the internal shocks in the jet. In addition, the generated/amplified magnetic field 
in the shock may survive for only a short time, this will further decease its value. We set 
/sh = 0.1 in our model. Fortunately this value is not very important since it can be absorbed 
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Fig. 1. — Spectral modeling results for XTE J1118+480. The fit was made with a model 
consisting of an inner hot accretion flow and an outer cool thin disk. The parameters of the 
model are r^r = 300rs, Mq = O.OSMEdd; « = 0.3, j3 = 0.9,6 = 0.5. The dashed line shows 
the emission from the inner hot accretion flow, the dot-dashed line shows the emission from 
the outer cool disk, and the solid line shows the sum of the two. The model explains the 
EUV and X-ray data quite well, slightly under-predicts the optical/UV, and significantly 
under-predicts the IR and radio fluxes (the radio data are shown in Fig. 2). Note that two 
sets of EUV data are shown, for two different choices of Nh. The X-ray spectral break at 
~ 10^^ '' Hz may not be real (see text for details). 
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Fig. 2. — Accretion-jet model of XTE J1118+480. The dashed and dot-dashed hnes show 
the emission from the hot and cool accretion flows, respectively, as in Fig. 1. The thin solid 
line shows the emission from the jet. The sum of the three components, shown by the thick 
solid hne, explains the spectrum all the way from radio to X-rays. The dotted line shows the 
synchrotron emission from power-law electrons that might be present in the hot accretion 
flow. 
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Fig. 3. — Advection factor / (defined in eq. 3) and the scaled mass accretion rate, M(r) / Mq, 
as a function of radius for the hot accretion flow model shown in Fig. 1. Negative values of 
/ indicate that the accretion flow is in the "LHAF" regime rather than the ADAF regime 
at these radii. The solution is basically an ADAF. 
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Fig. 4. — Model spectra from the hot accretion flow for three choices of Mq/M-^^^: 0.04 
(dotted), 0.05 (dashed), and 0.06 (dot-dashed). All other parameters are held fixed. 



